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Abstract 

The local Yang-Baxter equation (YBE), introduced by Maillet and 
Nijhoff, is a proper generalization to 3 dimensions of the zero curvature 
relation. Recently, Korepanov has constructed an infinite set of inte- 
grable 3-dimensional lattice models, and has related them to solutions 
to the local YBE. The simplest Korepanov's model is related to the 
star-triangle relation in the fsing model. In this paper the correspond- 
ing discrete equation is derived. In the continuous limit it leads to a 
differential 3d equation, which is symmetric with respect to all permu- 
tations of the three coordinates. A similar analysis of the star-triangle 
transformation in electric networks leads to the discrete bilinear equa- 
tion of Miwa, associated with the BKP hierarchy. 

Some related operator solutions to the tetrahedron equation are 
also constructed. 
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1 Introduction 



The Yang-Baxter equation (YBE) |T], |j plays the central role in the for- 
malism of the quantum inverse scattering method of solving integrable 2d 
quantum systems Q. The corresponding dynamical equations can be repre- 
sented as a zero-curvature condition (Lax equation) for some auxiliary linear 
problem. 

The notion of the local YBE has been introduced by Maillet and Nijhoff 
in ||], U as a proper 3-dimensional generalization of the zero-curvature con- 
dition. Korepanov in |], has pushed forward this idea by constructing an 
infinite set of discrete 3-dimensional integrable systems, and relating them 
to solutions to the local YBE, the simplest solution being equivalent to the 
star-triangle relation (STR) in the 2-dimensional Ising model. 

In this letter first we consider the star-triangle transformation in electric 
networks, and show that, being interpreted as the local YBE, it leads to the 
integrable discrete bilinear equation, first discovered by Miwa in ||. This 
is the content of Sect. Then, in Sect. [3|, in the same way we analyse 
the Korepanov's model associated with the STR in the Ising model. The 
corresponding integrable discrete equation in the continuous limit leads to 
the 3-dimensional non-linear differential equation, the symmetry group of 
which contains all permutations of three coordinates. In Sect. ^jsome opera- 
tor solutions to the tetrahedron equation, related with the abovementioned 
models, are constructed. 

2 The star-triangle transformation in electric net- 
works and Miwa's equation 

The star-triangle transformation in electric networks is historically the first 
transformation of this kind known since the last century. It is given by the 
following formulae: 

rjr'j = r[r 2 + r 2 r 3 + r 3 r[= r 1 r 2 r 3 /(r 1 +r 2 + r 3 ), j = 1,2,3, (2.1) 

where the resistances {r'j} and {rj} correspond to the star and triangle, 
respectively. There is a STR in the usual sense ||, underlying this transfor- 
mation. Define a Gaussian exponent: 

w(r\q) = exp(-iq 2 /2r), q,r E R. (2.2) 
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Then, it is easy to check that the following STR holds: 

J II I I \ ll\ I I \ ll\ I I \ II I 

dq w(r 3 \q-q )w{r 2 \q )w[r x \q -q 



J 27rr[r' 3 /ir2w(ri \q)w(r2 \q — q')w{r^\q'), (2-3) 



where the parameters {r'-}, {rj} satisfy relations (|2.1|) . The STR can be 
rewritten in the operator form as follows. 

Let C\p,q] be the Heisenberg algebra generated by Hermitian operators 
p and q, satisfying the Heisenberg permutation relation: 

qp — pq = i. (2-4) 



Then, interpreting the variables q,q',q" in the STR ( |2.3| ) as the eigenvalues 
of the operator q, one can write the STR in the form: 

w(—l/r 3 \p)w(r 2 \q)w( — l/r[\p) = w{ri\q)w( — l/r2\p)w{r 3 \q), (2-5) 

where the function w(r\q) is defined in d2.2|). A particular choice of the 



parameters, r[ri = —1, n = r 3 = —T2 = ±1, in (2J3) gives a constant 



solution to the YBE, which is discussed in [11| in the context of the quantum 
discrete sine-Gordon model. 

To interpret ( p.5|) as a local YBE, introduce a linearly independent set of 
vectors in it! 3 , {ei, e^i ^3}, and the 3-dimensional lattice, generated by these 
vectors with integer coefficients: 

M = {n = niei + n 2 e 2 + n 3 e 3 | ri\, n 2 , n 3 G Z} (2.6) 

Let ai,a 2 ,a3 be the functions on M: 

af M B n 1— > aj(n) G C, j = 1, 2, 3. 

The form of the identity ( p.5[ ) suggests the following parametrization: 

r[ = -l/ai(n), r' 2 = -a 2 {n - e x - e 3 ), r' 3 = -l/a 3 (n), 

ri = ai(n - e 2 - e 3 ), r 2 = l/a 2 (n), r 3 = a 3 (n - ei - e 2 ). 



Substituting these expressions into (|2.5D, we obtain the local YBE as is 
defined in g §: 

w(a 3 (n)|p)t(;(-a 2 (n - ei - e 3 )|<?)u>(ai(n)|p) 
= w(ai(n - e 2 - e 3 )\q)w(-a2(n)\p)w(a 3 (n -e\- e 2 )\q). (2.7) 
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The star-triangle transformation formulae ( |2.1| ) acquire now the form of 
difference equations: 

aj(n -e k - ej) __ 1 + a 2 {n -e\- e 3 )(ai(n) + a 3 (n)) , . , _ , 

7 \ — i \ ( \ i \3i fc i £ J — I 1 ' A "X- 

aj(n) ai(n)a 3 (n) 

(2.8) 



Theorem 1 T/ie substitution 

a ( n = a — — . . . — {j, fc, i} = {1, 2, 3}, 

J r(n + ei + e 2 + e 3 )r(n - e 3 -J 



where aj are some numerical coefficients, reduces system of equations (2.8) 
to the only equation on the function r(n): 

4 

y } j a.jT(n + ej)r(n — ej) =0, e± = —e\ — e 2 — e 3 , «4 = aia^o^. (2.9) 
i=i 

The proof is through the straightforward substitution. Miwa in || has 
shown that equation ( |2.9| ) in the suitable continuous limit reduces to the 
BKP equation, while it itself is equivalent to the whole BKP hierarchy. The 
limit e — ► 0, where ctj = (3je, j3j being fixed, reduces equation (|2.9|) to the 



Hirota's equation of the discrete Toda system [12]: 
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Y / P j T(n + e j )r(n-e J )=0. (2.10) 

3 A difference equation associated with the STR 
in the Ising model 

In this section we make a similar analysis to that of sect. ^ of another STR, 
associated with the Ising model. The corresponding dynamical system has 
been introduced and studied by Korepanov in || f7), where he has proved 
that this system has enough set of conserved quantities to be integrable, and 
described how to construct solutions of the finite gap type. 

Define the Bolzmann weights of the Ising model as a function W: 

W:C x {1, -1} 3 (x, a) ■-> W{x\o) = l + x + a-xaeC, (3.1) 
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The STR has the form of the following identity: 

]T W(x' 3 \aa")W(x' 2 \a")W(x' 1 \a"a') 
a"=±l 

= (l + x' 1 x' 2 x' 3 )W(x 1 \a)W(x 2 \aa')W(x 3 \a'), (3.2) 

where 

xjx k = (x'jx'f, + Xi)/(1 + x' x x' 2 x' 3 ), {j, k, 1} = {1, 2, 3}. (3.3) 



To rewrite STR ( |3.2j ) in a matrix form, first define the algebra C[Z, X] of 
2-by-2 matrices, generated by matrices Z and X, satisfying the relations: 



z 2 = X 2 = 1, zx = -xz. 



With these definitions, the STR (3.2) can be rewritten in the form (see also 
63): 

W((l - 4)/(l + x' 3 )\X)W(x' 2 \Z)W((l - x[)/(l + x[)\X) 

= RW( Xl \Z)W{(l - x 2 )/(l + x 2 )\X)W(x 3 \Z), (3.4) 



where 



R ^ (i + x 2 )(i + x , 1 44) 



(1 + ^(1 + 4) 

Next, introduce the new parameters through the formulae: 

4£j = (xj - 1/ Xj ) 2 , 4/$ = (x'j - l/x'rf, j = 1,2,3. (3.5) 

Consider three complex valued functions bj{n), j = 1,2,3, on the 3d lat- 
tice M defined in (|2.6|), and identify their values with the above defined 
parameters as follows: 

fi = l/6i(n), £2 = b 2 (n- e x - e 3 ), £ 3 = l/b 3 (n), 

ei = l/6i(n-e 2 -e 3 ), £ 2 = b 2 (n), = l/b 3 {n - e x - e 2 ), (3.6) 
Operator identity ( |3.4j ) now takes the form of the local YBE. 

Theorem 2 Under the substitution 

= jfe±affi±3) -,*, i} = {1,2,3}, 

/(n + ei + e 2 + e 3 )/(n - e,) 
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the star-triangle transformation formulae 3Q together with (3.5) and (3.6) 
imply the following equation on the function f{n): 



4 4 4 

2 £(/(n + e,)/(n - e,-)) 2 + 4 JJ /(n + e,) + 4 J] /(n - e,) 

V 

^/(n + e i )/(n-e i ) , e 4 = -ei - e 2 - e 3 . (3.7) 
V =1 / 

To prove the theorem one has to exclude the dependent variables by cal- 
culating the corresponding resultants, which in intermediate steps are cum- 
bersome. In these calculations the Maple computer system has been used. 

Consider now a particular continuous limit of the equation Q3.7j ), Let e 
be the lattice spacing. Introduce the field 



(era) = 31og(/(n)), n G M. (3. 



Theorem 3 In the limit e — > with x = (^1,^2,^3) = era € i? 3 fixed, the 
equation ({3. 1 ) for the function ( \3.8j ) is reduced to the following differential 
equation: 

6A 2 0][>f0) 2 = 2(A 3 0) 2 + (A 2 </>) 3 + 8]T(^) 3 , (3.9) 
i=i 3=1 

where <fi = 0(x), x = (xi, x 2 , X3), 

4 

A a = £c£, a = 2,3, 
i=i 

9,- = d/dxj, j = 1, 2, 3; <9 4 = -di - d 2 - d 3 . 

To prove the theorem one has to expand both sides of (|3.7| ) in power series 
of e up to the order of e 6 . 

Thus, we have obtained a differential equation which should be inte- 
grable as a continuous limit of the discrete integrable equation, and which is 
symmetric with respect to any permutation of the three coordinates in i? 3 . 
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4 Solutions to the tetrahedron equation 



The quantum theory of discrete systems, associated with the local YBE, 
should be described somehow by solutions to the tetrahedron equation (TE). 
Here we show how one can construct certain operator solutions to the TE 
using the star-triangle transformations considered in sections || and ^. 

First, consider the case associated with electric networks. Let i?+ be the 
set of positive real numbers: 

R + = { r e R | r > 0}, 

and i be the standard involution in it: 

l:R + ^R + , i*(r) = l/r, (4.1) 

where the variable r is interpreted as a coordinate function on R + . Define 
the following rational involution in R$: 

V :R 3 + ^R%, (p *(r j )=r j n+r2 + r3 , j = 1,2,3, (4.2) 

nr 2 r 3 

where r\ , r 2 , r$ are the coordinate functions of the corresponding multiples 
in R^. Note, that tp commutes with the action of the permutation group on 
R^_ . Formulas (^^) are derived in fact from the star-triangle transformation 
( |2.1|) through the following identification: 

r^l/^fo), j = 1,2,3. 

Define now one more involution of R+ as the following composition: 

<S> = L 2 oipaL 2 ;R%^R% y (4.3) 

where t 2 = id x l x id. Let $i 7 7 c , 1 < i < j < k < 6, be the involutions of 
R+, which act as on the subspace R+ C R\, specified by the i-th, j-th 
and k-ih -R+-multipliers, and identically, on the others. 

Theorem 4 The following TE holds in R^_ : 

$123 $145 $246 ° $356 = $356 ° $246 ° $145 ° $123- (4.4) 



The proof is the straightforward verification. The identity (fO|) implies the 
operator identity 

$123$i45$246$356 = $356$246$145$^3 (4-5) 
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in the linear space J-(R\) of complex- valued functions on R^_. Thus, we 
have obtained an operator solution to the TE. Solutions of such kind, called 
solutions to the "functional" TE, were considered also by Korepanov j!3| |. 

Let us comment on the limit, corresponding to Hirota's equation ( |2.10| ), 
For A S R + define the "rescaling" mapping: 

gx:Rl^Rl 



9\{ r j) 

and consider the following composition: 



Ar 



J = 1,2,3, 



$a = gi/\°$°gy-R + R+, 

where $ is defined in ([OP , flO|), ([O]). Clearly, $ A satisfies the TE (gT 
for any A 7^ 0. The limit A — * 0, however, is non-singular: limA_>o &\ = 
and the mapping $0, associated with the Hirota's equation (2.10), satisfies 
flOp as well. 

Consider now the star-triangle transformation fl3.3|) , corresponding to 
the Ising model. Let / be the open unit interval in R: 

I = { x G R I < x < 1}, 

and j be the following involution in I: 

f(x) = (l-x)/(l + x), (4.6) 

where x is considered as a coordinate function on /. Define the following 
involution in I 3 : 

V>:/ 3 -» I 3 , 



(1 - x k x{) 2 - x 2 Ax k - Xl) 2 



{j,k, 1} = {1, 2, 3}. (4.7) 



^ (1 + x k x{) 2 - x 2 (x k + Xlf 
This mapping is related with ( |3.3[ ) through the following identification: 

r(x j ) = (l-x' j )/(l + x' j ). 

The mapping 

^= j2 o^o j2 :/ 3 ^/ 3 , (4.8) 

should satisfy the TE identity (4.4) with $'s replaced by ^'s and considered 
in J 6 . We have verified this identity numerically, and claim it as a conjecture. 
The corresponding mapping 

^*:F(I 3 ) -» F(I 3 ) 
can be interpreted as a linear operator in the functional space J- (I 3 ). 
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5 Summary 



Integrable discrete equations in three dimensions can be obtained through 
solutions to the local YBE. Particularly, the star-triangle transformation in 
electric networks (2.1) can be associated with a solution to the local YBE 

I 



which leads to the integrable bilinear Miwa's equation (2.9). For the 
Korepanov's dynamical model, related with the star-triangle transformation 
in the Ising model, in this way we obtain another discrete equation (3.7). 
The continuous limit of the latter is described by a three dimensional differ- 
ential equation Q3.9] ), which is symmetric with respect to any permutation 
of coordinate axes. 

The operator solutions flOD , Q , Q , and ((D|, Q, flO| to the 
TE equation ( [4.4| ), connected with these dynamical models, may help to 
quantize the latter. Note, that the quantum version of Hirota's equation 
( |2.10| ) has been constructed in [14] in the context of the quantum inverse 
scattering method. 
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